Theory of a single-atom laser including light forces 



in 
o 
o 

Oh- 
< 
(N 



> 

o ■ 
^t- : 
o . 
in ■ 
o ■ 

^ : 
Ok 
^— > ■ 

G ■ 
cd 

3 : 
cr 



Thomas Salzburger 1 , Peter Domokos 2 , and Helmut Ritsch 1 
1 Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Technikerstrae 25, Austria 
2 Research Institute for Solid State Physics and Optics, P. 0. Box 49, H-1525 Budapest, Hungary 

We study a single incoherently pumped atom moving within an optical high-Q resonator in the 
strong coupling regime. Using a semiclassical description for the atom and field dynamics, we 
derive a closed system of differential equations to describe this coupled atom-field dynamics. For 
sufficiently strong pumping the system starts lasing when the atom gets close to a field antinode, and 
the associated light forces provide for self-trapping of the atom. For a cavity mode blue detuned 
with respect to the atomic transition frequency this is combined with cavity induced motional 
cooling allowing for long term steady-state operation of such a laser. The analytical results for 
temperature and field statistics agree well with our earlier predictions based on Quantum Monte 
Carlo simulations. We find sub-Doppler temperatures that decrease with gain and coupling strength 
and can even go beyond the limit of passive cavity cooling. Besides demonstrating the importance 
of light forces in single-atom lasers, this result also gives strong evidence to enhance laser cooling 
through stimulated emission in resonators. 

PACS numbers: 32.80.Pj, 42.50.Vk, 42.50.Lc 



I. INTRODUCTION 

A single incoherently pumped atom within a high-Q 
optical resonator constitutes the smallest and conceptu- 
ally simplest conceivable laser pj . On the one hand, this 
model is particularly interesting as it allows to theoreti- 
cally study important aspects of laser physics analytically 

HH, while > 

on the other hand, there is great technolog- 
ical interest in very small tailored coherent light sources. 
Already two decades ago laser like systems have been set 
up in the microwave regime 0, Q . With the tremendous 
recent progress in laser cooling and micro cavity technol- 
ogy such systems have now indeed been experimentally 
realized [aHSQ m the optical regime. This requires ul- 
tracold atoms trapped in a rather small volume between 
mirrors of extremely high quality. In this regime the light 
forces induced by the cavity field on the atom get impor- 
tant and have to be accounted for. In a first approxima- 
tion these forces are detrimental by heating the atomic 
motion and limiting the operation time of the system. 
This heating is less problematic for a trapped ion, but 
here the spatial requirements for the ion trap prevent to 
reach the strong coupling regime. In optical cavity QED 
setups such heating is significant and strongly shortens 
the interaction time |10| . 

In a recent paper we have shown that for carefully cho- 
sen operating conditions one can reverse the detrimental 
effect of heating and combine gain with optical cooling 
and trapping yjj. In this way the laser field generated 
by the atom can be used to simultaneously trap and cool 
the atom, which leads to a self-sustained sustained laser 
operation. Interestingly, the temperature attained by the 
atom can be even lower than for free space Doppler cool- 
ing at comparable parameter values. This stems from 
the strong nonlinear dependence of the field intensity on 
the atomic position, which gives new prospects of devel- 
oping a novel laser cooling method enhanced by stimu- 
lated emission. The numerical results obtained in Ref. Ill 



clearly demonstrate the significance and potential useful- 
ness of light forces for single-atom lasing [lfj, but they 
only provide little insight in the basic physics going on 
to facilitate this interesting behavior. 

In this work we develop a systematic semiclassical de- 
scription of such a single-atom laser in the spirit of the 
successful models developed for cavity cooling Ap- 
proximating the atomic center of mass motion by a point 
particle, we are able to find a set of coupled equations for 
the internal atomic dynamics and field evolution param- 
eterized by the atomic position. In contrast to laser cool- 
ing we cannot adiabatically eliminate the upper atomic 
state as we need an inverted atom for gain. Neverthe- 
less, by use of a factorization approximation for higher 
order atom-field expectation values, we still can derive a 
closed system of coupled equations for the combined dy- 
namics. As a proven but simple model for pumping we 
use an inverted heat bath approach developed by Haken 
several decades ago [12j. From these internal dynamical 
equations we are then able to derive the average light 
potential as well as friction and momentum diffusion co- 
efficients for the atomic center of mass motion in an ana- 
lytic form. This allows detailed studies of atomic motion 
and laser field evolution that can be checked and com- 
pared to quantum Monte Carlo simulations for selected 
parameters. 



II. MODEL 



Let us consider a single inverted two-level atom freely 
moving in the field of an optical resonator with high 
finesse. The relevant mode with frequency lu c = kc 
is detuned from the atomic transition frequency u a by 
A = u> c — u> a . The dipole-resonator-field interaction 
Hamiltonian in the rotating wave and the electric dipole 
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approximation (H = 1) then reads 

H = -Aer + er_ + iG (aV- - cr+a) (1) 

which is written in a frame rotating at the cavity res- 
onance frequency. Here <r_ (o + ) and a denote, re- 
spectively, the atomic lowering (raising) operator and 
the bosonic field operator for the cavity mode while 
the position dependent atom-field coupling is given by 
G = gcos(kx). 

Both the field and the atom are coupled to the en- 
vironment, which is modeled by Markovian decay pro- 
cesses with rates 2k (photon loss via the mirrors) and 
27 (spontaneous emission). Using standard techniques 
of quantum optics, we can derive the following master 
equation: 

^£?= -i[H,g] +C K g + C 1 g. (2) 

It describes the time evolution of the resulting open sys- 
tem including decay of the resonator mode (C k q) and the 
atomic upper state (C 7 g). 

In order to feed energy into the system, the atom is 
driven externally by incoherent excitation at rate 2v. As 
a simple but still quantum mechanically consistent way to 
incorporate such a pumping mechanism, we model it by 
inverse spontaneous emission. This proven method has 
been introduced already in the early quantum models of 
lasing and we simply have to add a corresponding 
Liouvillian term to Eq. 

C V Q = v(2o+Q<J- — (J-CF+Q — QO-O+) . (3) 

In fact it is largely equivalent to pumping from the 
ground state to a third intermediate level, with a fast 
incoherent decay to the upper atomic state. 

The above master equation is equivalent to the follow- 
ing set of Heisenberg-Langevin equations 13] : 

d = —na + G<7_ + t; a (4a) 
a- = (iA - 7 - v)o- + Ga z a + £ CT (4b) 

& z = -2(7 + v)a z - 2G(a t cr_ + a+a) + 2{v - 7) + ■ 

(4c) 

Here we have introduced noise operators originating from 
the coupling of the system to the environment. While 
their expecation values vanish when the environment is a 
T = heat bath, their nonvanishing correlation functions 
are given by: 

-t)> = 2««(t) (5a) 
(Utn(t-r))=2 1 5(r) (5b) 
-r))=2i/a(r) (5c) 
(6(*)&(t - t)) = (4( 7 - u)(a t ) + 4( 7 + u))S(t) . (5d) 

In this paper we assume that the kinetic atomic tem- 
perature stays well above the recoil limit ksT rt , c = 



H 2 k 2 /(2m), which of course has to be self-consistently 
checked at the end. This allows for a so called semiclas- 
sical approximation, where the particle's position and 
momentum are treated classically and enter the equa- 
tions for the field and internal atomic dynamics simply 
as real parameters x and p. A systematic way to derive 
these equations for an atom in a cavity field is, e.g., pre- 
sented in Ref. |llj. In this approximation simulations of 
atomic trajectories will then be governed by the following 
Langevin-type equations 

x = p/m (6a) 
P = F + £ , (6b) 

where F denotes the average force acting on the atom 
and £ is a noise term giving rise to momentum diffusion. 
Both of these values have to be calculated from the cor- 
responding solution of the master equation Eq. J5J • 

III. LAMB SEMICLASSICAL MODEL 

Let us first try to get some qualitative insight into the 
dynamics of our system and neglect all the noise terms 



£i and replace all operators by c-numbers, i.e: (a) = a, 
(<r_) = s, and (o z ) — z. Hence, we get the following set 
of coupled differential equations: 

a = —na + Gs (7a) 
s = (iA - 7 - v)s + Gza (7b) 
i = -2(7 + v)z - 2G(a*s + s*a) + 2(y - 7) . (7c) 



This can be readily solved for the steady state for an 
atom at fixed position. Although any operator and noise 
correlations are neglected in this "Lamb-type" model, we 
still get some vital insights into the dynamical character- 
istics of our system. As an immediate consequence, Eqs. 
l |7| l imply the continuity equation 

v{l-P)= 1 P + nN (8) 

which describes the energy balance in the system due 
to pumping and losses via the atom and the cavity in 
the stationary state. Here we used the atomic ground 
and excited state populations 1 — P = (1 — z)/2 and 
P = (1 + z)/2 and the photon number N = \a\ 2 . This 
is a universal relation independent of the particle's posi- 
tion and will be recovered several times throughout the 
paper. Linking P to N, Eq. (JSJ immediately yields the 
atomic population from the intracavity intensity which 
will appreciably simplify the analysis in the following. 

Obviously for an atom fixed at a node of the cavity 
field where the atom-field coupling strength vanishes, 
the photon number is zero as well and the atomic up- 
per state population is P = vj(y + 7). When the atom 
moves into regions where G exceeds the threshold value 
G t h — v^CCr + v ) 2 + A 2 )/(^ 7): the atom-field cou- 
pling opens an additional decay channel via the cavity 
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mode. Indeed, one can calculate the rate of emission 
into the cavity mode 



W 



(7 + ^)G 2 
(7 + vY + A 2 



(9) 



The general behavior of the cavity photon number is 
depicted in Fig.^ a where we have plotted N (solid line) 
as a function of the atomic position along the cavity axis 
within half a wavelength. Note that the cavity field starts 
to be populated with G crossing a threshold value in a 
highly nonlinear way. According to Eq. JSJ), this sudden 
increase has to be accompanied by a corresponding drop 
in the atomic population inversion z — k/W (dashed 
line). At the same time there will be a big change of the 
light force on the atom, originating from the modified 
optical potential. If the atom is a high field seeker, this 
already lets one expect a possible tight confinement of 
the atom. 

For an atom at rest at a fixed position the mean force in 
steady state is simply proportional to the photon number 
as well as the gradient of the mode function and explicitly 
reads: 



F 



2kA VG 

j + v G 



N . 



(10) 



Obviously, F will be zero at antinodes, where N is maxi- 
mal, due to the vanishing gradient of the mode function. 
Notice that for G — ► 0, N tends faster to such that ex- 
pression i|10[) remains welll defined zero. Since W, z, and 
N are even functions of the detuning A above thresh- 
old, we get F(— A) = —F(A), and the atom will be a 
high-field seeker for A > 0. 

This can be seen in Fig.^b which shows F (solid line) 
as a function of x for A = 200k as well as the correspond- 
ing light potential, both in arbitrary units. Although this 
sounds contradictory first when compared to standard 
formulas for the optical potential, one has to remember 
that the atom is inverted and thus the sign of the light 
potential is reversed and dominated by the upper level 
Stark shift. Note that within this approximation a fixed 
atom will not feel any mean force unless G > Gth- 

Let us now look at the full coupled dynamics of atom 
and field by simultaneously integrating Eqs. © and (7J), 
which can be easily performed numerically. In Fig. ^ 
we show the evolution of the particle's position x and 
momentum p as well as the intracavity photon number ./V 
for a typical set of parameters where one gets trapping. 
The atom starts at some random position (plotted in 
units of A) initially moving fast along the cavity axis. 
Gradually, its motion gets damped until its kinetic energy 
falls below the potential depth, and the atom is then 
confined to oscillate in a single well. Here its kinetic 
energy is still reduced further but at a much slower rate. 
As a remarkable feature the photon number suddenly 
undergoes a drastic increase as the atom gets trapped. 
This is due to the fact that the atom remains close to 
antinodes and never enters a spatial region, where the 



laser threshold is not fulfilled. This could experimentally 
clearly be used to observe trapping in real time. 

The physical mechanism responsible for the fast dissi- 
pation of the particle's motional energy during the initial 
stage is yet another variant of Sisyphus cooling. When- 
ever the atom enters a spatial region where the system 
falls below threshold, the momentary photon number is 
higher than at the time later when it reenters the lasing 
region. Hence, it is pulled back stronger during leaving 
than sucked in during reentering, which gives rise to net 
friction forces. This type of friction force even continues 
while the particle is oscillating in a single well as illus- 
trated in Fig. ^ d- It shows a cutout of the trapping 
phase. Due to the finite response time of the cavity field 
to the atomic position the cavity field attains a maximum 
value every time shortly after the atom has passed a field 
antinode (see 1) and the atom still feels a significant fric- 
tion. On the other hand, when the atom approaches a 
turning point (2), the intensity reaches a minimum re- 
sulting in a smaller accelerating force towards the center. 
This is of course very similar to passive cavity cooling in 
principle. However, here the response time of the field is 
not only dominated by the cavity decay rate but by the 
full laser field dynamics. This can strongly enhance the 
effective friction and lower cooling time and the steady- 
state temperature. 



IV. QUANTUM RATE EQUATIONS FOR THE 
INTERNAL DYNAMICS 



Let us now go beyond the simple factorized c-number 
model and include fluctuations due to the interaction of 
the system with the environment represented by the noise 
operators in Eqs. J3J. As one central consequence, these 
operators induce momentum diffusion of the atomic mo- 
tion counteracting the cooling process and prevent the 
atom from stopping completely at a field antinode. They 
also introduce fluctuations in the photon number and 
atomic occupation probabilities. 

In the following quantum model we include the corre- 
sponding noise terms but we still assume a rather local- 
ized atomic wavepacket. This allows to replace atomic 
momentum and position operators by their average val- 
ues, i.e. we treat the atom like a classical Brownian par- 
ticle in the optical potential. However, we will keep the 
quantum correlations between the cavity field and the 
atomic polarization which were neglected in the previ- 
ous section due to factorization. We therefore base our 
treatment on second order operator products which turn 
out to obey a closed set of equations. Using the abbre- 
viations <f> = a^a, n = er_|_er_, X = crcr_ + c+a, and 
A = {a)a- — <j+a)/i we get: 
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FIG. 1: a Photon number JV (solid line) and population in- 
version 2 (dashed line) for the steady state as a function of 
the particle's position x. The parameters are (7, v,g,A) = 
(10, 20, 100, 200)k. 

b Stationary force acting on the atom F (solid line, a. u.) and 
corresponding potential U (dashed line, a. u.) for the same 
parameters. 

c Time evolution of the photon number N, particle position 

x (units of A = 2n/k), and particle momentum p (normalized 

to be unity when the atomic temperature equals the Doppler 

temperature T — h-y ) for the same parameters. 

d Cutout of c demonstrating the cooling mechanism (see 

text). 



$ = - 2/t$ + G£ + T<e, (11a) 
II = - 27II - GE + 2z/(l - II) + T n (lib) 

s = - rs - aa + ?;g[s, a] + r s (lie) 

A = -rA + AS + T A . (lid) 

Here r = K + 7 + ^is the total damping rate of £ and A 
that gives the atom field interaction energy and is closely 
related to the force. The quantum fluctuations are con- 
tained in the operators Tj which again are fully deter- 
mined by their second-order correlation functions given 
in appendix [S] Note that Eqs. (|ll(l are exact but still 
clearly nonlinear. Thus, the corresponding equations for 
their expectation values are not closed and we have no 
explicit solution for their steady state. The difficulties 
arise from the operator product in Eq. (|llc|) . An extra 
equation for ([£, A]) of course will inevitably incorporate 
higher order operator products resulting in an infinite 
hirarchy of equations. Following an idea developed in 



earlier laser models 24] , we break this loop by replacing 

i[E,A] = 2(2n-l)$ + 2n«2Z$ + 2n, (12) 

where Z = (211 — 1) is a real parameter that later can be 
calculated from Eq. |JSJ self-consistently. This approx- 
imation means that we drop part of the quantum cor- 
relations between the atomic populations and the field 
intensity. Fortunately, this turns out to play a minor 
role in the calculation of the system- variable expectation 
values in the parameter regime we are interested in. The 
factorized equations for S = (<&, II, £, A) now read 



— 3 = M3 

dt 



(13) 



where we have defined v = (0, 2v, 0, 0) and 
/ 

M = 



-2k 
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-2(7 + f) 


-G 
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2G 


-r 


-A 
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-r 



(14) 



If the coupling strength g is much less than the damping 
rate T or the detuning A, the operators £ and A will 
adiabatically follow the values of <& and II such that one 
is allowed to eliminate them. Adiabatic elimination of 
£ and A then yields the quantum rate equations for the 
photon number N — ($) and the particle's upper state 
population P = (II), 



N = -2(« - ZW)N + 2WP 

P = -2(7 + W + v)P - 2ZWN - 



2v . 



Again, 



W = 



TG 



r 2 



(15a) 
(15b) 



(16) 



denotes the emission rate into the resonator. In contrast 
to Eq. W now comprises the combined rate T. It 
tends towards the rate found in the previous section when 
the cavity relaxation time k~ 1 is much longer than any 
other timescale of the system. 

Apart from relation JSJ that immediately follows from 
the rate equations in steady state, there are two main ef- 
fects arising from the strong atom-field coupling. (Note 
that © is recovered whether or not the factorization 
(|12fl is made.) First, the emission of photons into the 
resonator field is represented by the terms explicitely in- 
terlinking Eqs. H15(l . While the atomic population P de- 
creases, the field intensity grows due to the source term 
2WP. 

Interestingly, both the atomic and the cavity linewidth 
are effectively modified as can be seen from the remain- 
ing terms proportional to W. This intricately affects the 
decay properties of the whole system and, e.g. , reduces 
the resonator linewidth for an inverted atom. As a re- 
sult, even the stationary solution depends on the atomic 
position in a highly nonlinear way. 
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to be reduced for several gain atoms in the mode where 
the inversion can be shared among many atoms. 
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FIG. 2: Up: stationary photon number obtained from the 
rate equations as a function of the atomic position x and the 
detuning A. The parameters are (j,i/,g) = (20, 25, 20)k. 
Center: atomic upper state population for the same parame- 
ters. 

Down: The dipole force F (plotted in arbitrary units) is an 
odd function of A. 



Let us emphasize here that the cooling limit for con- 
ventional cavity cooling is usually related to the cav- 
ity linewidth k 0]. An effectively gain reduced cav- 
ity damping rate thus can give rise to even lower final 
temperatures |15) . On the other hand, for an inverted 
atom momentum diffusion due to spontaneous emission 
is strongly pronounced and since W can be of the order of 
7, heating through dipole fluctuations will tend to raise 
the particle kinetic energy. This effect can be expected 



V. FORCES 

In the following we will investigate the atomic motion 
in more detail. While the atom moves under the influence 
of the light forces induced by the cavity field, the atom 
modifies the light field dynamics according to its position. 
For a passive resonator this mutual influence is a well- 
known feature of c avity QED and has been already seen 
experimentally p"?l UM There arc further complications 
in the present system. First, the atom itself generates the 
light field it interacts with. Second, the photon creation 
from the incoherent pump is a highly nonlinear process 
exhibiting threshold. In particular there is no light in the 
mode without an atom close to an antinode. As the atom 
also provides gain inside the resonator, the lasing-induced 
trapping and cooling effects will be strongly enhanced 
when the system operates above threshold. 



A. Photons and forces in the steady state 

As is well known, the radiation pressure force on an 
atom at rest in a standing wave field cancels on aver- 
age and only the dipole force is left ,2(|. This remains 
true for a standing wave cavity field, and the only contri- 
bution to the net average force arises from the reactive 
response of the atom to the field dynamics (dipole force 
F). For a slow enough atom the state of the field dy- 
namically adjusts to its current position x and can be 
well approximated by the steady state for an atom fixed 
at x. In this adiabatic approximation we can calculate F 
from the stationary expectation value of the force oper- 
ator T that is given by the Heisenberg equation for the 
momentum operator V, 



in 



-VH = (VG)A. 



Inverting the matrix M and defining D 
get from Eq. lpHI)l : 



F 



2kvAG 
D 



(VG) 



(17) 
detM/4, we 



(18) 



Here the environment is considered at zero temperature 
and hence the noise operators have vanishing expectation 
values and do not contribute. From Eq. I|18|) we can read 
off that F is an odd function of the detuning A. By help 
of the solution 



N 



vYG 
D 



(19) 



for the photon number, we then see that Eq. i(T%|l tends 
to Eq. i|l(J|) for k <C 7 + v. This corresponds to resonator 
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fields close to coherent states with large mean intensity 
allowing a classical description. 

In Fig. we have plotted the stationary values of N, 
P, and F, respectively. As expected from Eq. the 
photon number is a monotonic function of the atom-field 
coupling G and decreases for growing mismatch between 
the atomic and resonator frequencies. In contrast to 
the Lamb model, even a slight displacement of the atom 
from a field node causes the atom to radiate into the las- 
ing mode although the system has not crossed the laser 
threshold. This small but nonzero mode occupation leads 
to a force acting on the particle now. Depending on the 
detuning, the atom is either pushed back to the node 
(A < 0) or attracted to the interaction region (A > 0). 
The resulting light potential shows minima at antinodes, 
where we would like the atom to be trapped for positive 
detunings, and hence we will concentrate on this situa- 
tion in the following. 

Fig-El a shows the potential depth V in units of the in- 
teraction energy as a function of the detuning A. In each 
curve the pumping rate v is adjusted in order to achieve 
constant maximum photon number N as indicated. The 
parameters are (j,g) = (10, 50)k. It is clear that V suc- 
cessively increases with growing photon number. On the 
other hand, large light intensities in the resonator will 
result in enhanced cooling and particle localization as we 
will see later. 



B. Friction and diffusion 

For finite particle velocities the field cannot follow the 
atomic motion instantaneously but will show some time- 
delayed response. Let us now derive the linear velocity 
dependence of the force, i.e. the friction coefficient (3. If 
the atom moves much less than a wavelength before the 
internal variables attain their stationary values accord- 
ing to a displacement of the atom, one can derive the 
linearized correction term (3 to the force as follows [l9j . 
Expanding the system operators in terms of the parti- 
cle velocity v and replacing the total time derivative by 
d/dt + wV we get a set of dynamical equations that can 
be solved systematically in different orders of v. Writing 
(H) w X° 



uX , the zeroth- and first-order variables 



obey 



X° = -M _1 v 

X 1 = -M _1 V fM" 



(20a) 
(20b) 



where the first line is precisely the afore discussed adia- 
batic solution. The friction coefficient is then given by 
(3 = (VG)(A 1 ) which is a rather lengthy expression and 
can be found in appendix [51 Note that in principle it is 
straight forward to derive the solution to higher orders in 
v, which would certainly produce more accurate results. 
However, here we are mainly interested in the parameter 
regime with a very low final temperature of the atom. In 
this limit the friction can be well approximated by the 
position- averaged force term linear in v. 
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FIG. 3: a Optical potential depth in units of the interaction 
energy for (7, g) = (10, 50)k. In each curve the pump rate v 
is adjusted resulting in constant maximum photon number N 
as indicated. 

b Position averaged friction coefficient in arbitrary units for 
the same parameters as in a. 

c The atomic equilibrium temperature can be lower than the 
Doppler limit To = ^7- 

d Ratio of atomic kinetic and potential energy. For JV > 1 it 
drops well below one when the cavity is far detuned from the 
atom indicating strong particle localization. 



In order to give an estimate for the temperature, we 
also need the momentum diffusion coefficient T> so that 
we can apply the Einstein relation poj 



k B T = (3/V. 



(21) 



Hence, we also have to calculate the force fluctuations 
due to the coupling of the system to the vacuum modes. 
In our case the cavity field as well as the atomic variables 
experience fluctuations around their stationary values 
which directly relate to force fluctuations. This counter- 
acts the cooling process and prevents the atom from stop- 
ping completely at antinodes. Eq. (|21|1 then describes a 
situation where the contributions of friction and heating 
cancel and the atom reaches an equilibrium momentum 
distribution. 

Similar to Brownian motion we calculate the diffusion 
coefficient T> from the linear growth term of the momen- 
tum spread due to field fluctuations [201. Here we use 
the approach first outlined in Ref. [21J which allows to 
approximately read off T> from the two-time covariancc 
of the force operator, 

{P{t)f{t t)) (Ht))(Ht r)) = 2V5(t) . (22) 

A more detailed derivation of T> can be found in ap- 
pendix [21 Here we only want to note that, in addition to 
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FIG. 4: a Atomic temperature as a function of the maximum 
photon number N in the cavity for different values of the 
pumping rate v. The other parameters are (7, g) = (10,50)fi; 
while the atom-field detuning A was continuously changed, 
b Corresponding ratio E/V . 



fluctuations of the cavity field, we have to consider mo- 
mentum diffusion owning to the random recoil of sponta- 
neously emitted photons and therefore have to add P rcc 
that can be found elsewhere [2(|. Note that in principle 
force fluctuations from the pumping mechanism would 
enter here as well. As we do not explicitly specify the 
corresponding mechanism and assume that the pumping 
occurs transversally, we will neglect this contribution at 
this point. 

Let us now come to some numerical examples. In Fig. 
[3]b we have plotted the position-averaged friction coef- 
ficient for the same parameters as in a. For very low 
intracavity fields we have /3 ra and the atomic motion 
is slightly accelerated rather than damped. Only above 
threshold a strong friction force arises. We see that for a 
given photon number we can have both heating and cool- 
ing, which already indicates that the system is active. In 
the regime where v < 7 the atomic population is not yet 
inverted and (3 > 0. For higher pumping strengths (i.e. 
larger A in Fig.[3J) there appears a change in the signs of 
the population inverstion as well as of the friction coeffi- 
cient. It turns out that Z > is, together with A > 0, 
a main condition to achieve cooling and, moreover, large 
atomic upper state populations imply low temperatures 
as we will see at the end of this section. 

Similar to the light potential also \(3\ shows a nonlinear- 
increase with the photon number. We can therefore ex- 
pect the cooling efficiency to be strongly enhanced for 
higher laser intensities. Indeed, the mean kinetic en- 
ergy of the atom continuously decreases when more and 
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FIG. 5: Atomic temperature as a function of the atom-field 
coupling constant g for different values of the pumping rate 
where (7, A) = (5,250)k. The marks show the results ob- 
tained from Monte Carlo wavefunction simulations. 

more photons are present in the resonator mode. This 
is demonstrated in Fig. [3] c where we have plotted the 
atomic steady-state temperature in units of the Dopplcr 
temperature Td = /17. We clearly get sub-Doppler cool- 
ing in the cavity field, which definitely proves the impor- 
tant role of the cavity for cooling and is the prerequisite 
to combine trapping and cooling. 

So far we have seen that the atomic motion can be ef- 
ficiently cooled when the atom is able to scatter ample 
photons into the resonator. An essential issue to achieve 
long term operation of such a device is strong particle 
localization at antinodes. Fig. [3]d shows the ratio of the 
atomic kinetic energy E and the optical potential V. Be- 
low threshold, this ratio is much larger than unity and 
the particle's position is almost evenly distributed along 
the cavity axis. For higher photon numbers, E/V can 
drop well below unity corresponding to strong localiza- 
tion. This is in big contrast to free space Doppler cooling, 
where the average kinetic energy is shown to be always 
larger than the optical potential depth. 

Like in most other cavity cooling schemes, the pumping 
strength has great influence on the trapping and cooling 
rates. Large light intensities induce fast and strong local- 
ization while the final temperature remains mostly unaf- 
fected there. However, in our system where the atom acts 
like a gain medium inside the resonator, also the parti- 
cle's kinetic energy shows strong dependence on the intra- 
cavity intensity. This is demonstrated in Fig.0]depicting 
the particle temperature and the ratio E/V vs maximum 
photon number N for different values of v. The large cor- 
relation of the internal dynamics and the atomic motion 
leads to situations where the atom is glued to antinodes, 
thereby radiating light into the resonator mode, which 
in turn carries away energy and entropy from the sys- 
tem via the cavity mirrors and thus decreases the atomic 
temperature. 

From Figs. |21 and 2| we see that our equations predict 
lasing, cooling, and trapping simultaneously. This occurs 
particularly for high photon numbers in the resonator 
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a v/A 

FIG. 6: a Atomic temperature in units of hn. For y > 1/2 
the temperature can fall below one which corresponds to the 
limit of passive cavity cooling. 

b Minimum temperature as a function of y = v/A and the 
parameter a determining the operating point of the system. 

and when the atom is far red detuned from the lasing 
mode. In this parameter range, though, the atomic up- 
per state has to be strongly populated and a large pump- 
ing rate v is required, which appears to be the central 
experimental bottleneck in this system. Naturally, this 
suggests to simultaneously use two or more atoms for 
gain. In this way they can collectively emit into the las- 
ing mode resulting in enlarged mode occupation without 
such stringent pumping requirements. We thus expect 
not only the threshold appear at lower pumping strengths 
but also advanced cooling and trapping. 

C. Comparison with numerical Monte Carlo 
simulations 

Since the early days of quantum optics, much atten- 
tion has been paid on the development of laser theories 
at different levels of sophistication. While most of the 
analytic models in the field of cavity- QED were based on 
rate equations at the beginning [23 , there has also been 
some work pointing out the shortcoming of the factoriza- 
tion approximation. More accurate models were intro- 
duced |23l |24| that could be applied in our case as well. 
We will, however, use a different approach here and di- 
rectly use numerical methods to solve our original master 
equation without approximations. Some first results on 
this were already published previously . Here we show 
a comparison with our analytical calculations in Fig. [S] 
Obviously, in the regime where lasing together with cool- 
ing coexists we find a surprisingly good agreement of the 
steady-state temperature and field expectation values. 

D. Good-cavity limit 

In an active system the cavity field response time is 
no longer simply given by the cavity decay k rate but 
gets dynamically modified. As k gives a lower limit on 
the kinetic temperature of the atom for passive cavity 



cooling, one could also expect changes here [l5j. I n the 
following we will study this in more detail and calculate 
the equilibrium temperature in the limit of very small k 
(particularly k <C v) by expanding friction and diffusion 
to first order in the decay rate k. In order to keep oper- 
ating conditions comparable, we scale the rate of photon 
emission into the resonator also linear in k, i.e. 

W = an. (23) 

Here the parameter a determines the operating point of 
the laser. Note the threshold condition k = W found 
in the classical model; hence a > 1 corresponds to the 
laser working above threshold. To keep the final expres- 
sion simple we further assume 7 ~ 0. This leads to the 
following rather simple expression for the atomic equilib- 
rium temperature 

t 2a 2 + (a - 1)V , s 

k B T = tin y - UL (24) 

2ay 

which depends on a and the amount of pumping ratio y = 
vj A. In Fig. El a we have plotted hsT/Un for different 
values of y. We find that for y > 1/2 it can drop below 
one and thus below the limit of passive cavity cooling. 
The respective minimum temperature, 

k B T = hK (vV + 2 - yj , (25) 

a monotonic decreasing function of y, is displayed in Fig. 
El b (solid line) . Again large pumping rates and hence 
large values of Z result in low temperatures. In addi- 
tion we show the corresponding a-parameter that remains 
slightly below one. Therefore the system operates a bit 
below threshold where the atom is mainly in the excited 
state and the interaction energy is very large. 



VI. CONCLUSIONS 

We presented a simple self-consistent analytical model 
for the coupled dynamics of an inverted atom described 
as a point particle moving in the field of an single-mode 
resonator. In good agreement with previously obtained 
predictions from Monte Carlo simulations we find that 
lasing, trapping, and cooling can simultaneously occur 
in such a setup, when the light mode is blue detuned 
from the atomic transition frequency and the pumping 
is sufficiently strong. This surprising result turns out to 
be closely related to the fact that an inverted atom is 
a high field seeker for blue detuning. Blue detuning is 
also a necessary condition for cooling as the missing pho- 
ton energy in the stimulated emission process has to be 
taken from the atomic kinetic energy. Luckily, this dissi- 
pation of kinetic energy via the resonator mode results in 
atomic equilibrium temperatures well below the Doppler 
limit and overcompensates the extra heating from the 
increased spontaneous emission of an inverted atom. 
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As a consequence high photon numbers not only imply 
stronger localization of the atom but also lower temper- 
atures. As an extra bonus the atom as a gain medium 
effectively can reduce the resonator field linewidth below 
the cavity linewidth, so that under favorable conditions 
temperatures even below the limit of conventional cav- 
ity cooling (ksT = Kk) are possible. This effect should 
definitely get more prominent for a larger atom number 
in the cavity. Hence, even for larger samples stimulated 
cooling could be connected with lasing in a combined 
atom-laser-photon-laser setup providing for a coherent 
atomic beam and light source. 

Let us finally remark that the fact that lasing is not 
necessarily connected to heating of the active medium 
but rather involves cooling could also prove important 



in rather different micro laser setups, e.g. on microchips, 
where heat production is a major issue preventing future 
miniaturization. 
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APPENDIX A: CORRELATION FUNCTIONS OF THE NOISE OPERATORS T; 

Analogous to the noise terms £j in Eqs. I@J, the operators contain only free input field operators and their 
expectation values vanish when evaluated at zero temperature. The remaining non-zero correlation functions are 

(T$(t)T$(t-r)) = 2kN5(t) 
(Tn(t)T n (t - r)) - (2 7 P + 2v{\ - P)) 5(t) 
(T E (i)T E (f - r)) = {2kP + 2 1 N + 2v{\ + N)) 6(r) 
(T A (t)T A (t-r)> = <T E (t)T E (t-r)) 
(T*(t)T n (t - t)) = (T n (t)T#(i - r)> = 
(T*(t)T E (t - t) + T E (f)T*(t - r)> = 2/s(E> S(t) 
(T»(t)T A (t - t) + T A (f)T*(t - r)) = 2k(A) 6(t) 
<T n (t)T s (t - r) + T E (t)T n (* - t)> = 2( 7 - i/)(E> <J(r) 
(TnWT A (t - t) + T A (t)T n (* - r)) = 2( 7 - j/)(A) 5(r) 
(T E (t)T A (< - r) + T A (t)T E (< - t)> = . 



APPENDIX B: FRICTION COEFFICIENT 

From the solution of Eqs. I|2(J|) we obtain the somewhat unhandy expression 



- G 3 r(4 K 2 ( 7 + vf T + G 2 ( 7 + v - k) (k 2 + (7 + ^) 2 ^)) (VZ) 

+ 2 K { (r 2 + a 2 ) (g 2 (« 3 - ( 7 + ^) 3 ^) - 2k 2 (7 + r) 

+ TG 4 (k - (7 + //)^) 2 + 2k (7 + v) Y 2 G 2 (k - (7 + } (VG) 



(Bf) 



APPENDIX C: DIFFUSION COEFFICIENT 



In the following we give a brief description of the calculation of the diffusion coefficient. Writing T = {£) + T, 
definition yields 

2D S(t) = (VG) 2 (T(t)T(i - r)) (CI) 
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since (T(i)) = 0. Here we see that the momentum spread directly arises from the noise exhibited by the interaction 
of the system with the environment via the operator T that, for quasi-stationary conditions, is given by 



T = — ((7 + v)AGZ T$ + kAG T n + /t( 7 + v)A T s + «( 7 + u)T T A + G 2 {k - (7 + v)Z) T A 
Inserting the correlation functions listed in appendix we find 



V = 



vG 2 {WGf 



£> 3 



2k 2 A 2 (7 - v + (7 + v)Z) (2k(7 + i/)r + G 2 (k - (7 + v)Z)) 
+ r 2 ( K 2 ( 7 + vf A 2 + [«( 7 + V )r + G 2 (« - ( 7 + v)Z)] 2 ) (1 + ± z) 
+ kTA 2 G 2 (2 K7 (A - Z) + K 2 + (7 + i/) 2 ^ 2 ) 



(C2) 



(C3) 



r 
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